function y = initial_function3(x)
y = 0;
for i = 1:30
    y = y + sin(i*x)/i*(1-(i/30)^2);
end